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A new method of quantum state tomography for quantum information processing is described. 
The method based on two-dimensional Fourier transform technique involves detection of all the off- 
diagonal elements of the density matrix in a two-dimensional experiment. All the diagonal elements 
are detected in another one-dimensional experiment. The method is efficient and applicable to a wide 
range of spin systems. The proposed method is explained using a 2 qubit system and demonstrated 
by tomographing arbitrary complex density matrices of 2 and 4 qubit systems using simulations. 

I. INTRODUCTION 



^ Quantum Computation offers exciting possibilities of solving complex computational problems 
^ using algorithms which exploit the quantum nature of the system. The idea, first proposed by 
iFeynman [1], is being feverishly pursued by many [2-6]. Several algorithms like Shor's factorization 
C3 algorithm, Grover's search algorithm, Deutsch-Jozsa algorithm, quantum Fourier transform, quan- 
tum counting and quantum error-correction codes have been developed and have clearly established 
• !the premise [7-24]. The last step in quantum information processing and quantum simulations is 
"the measurement of the output quantum state, known as quantum state tomography. In the case 
of ensemble systems this amounts to measuring the output density matrix. The output state of a 
quantum algorithm normally corresponds to some classical information and therefore, it is sufficient 
to measure all the diagonal elements of the density matrix which corresponds to the probabilities 
of various eigenstates. However, full quantum state tomography is generally carried out wherever 
possible, because of the following reasons; (i) knowledge of the full output density matrix allows one 
to find out the experimental errors and to calculate the fidelity of the implementation [5,10-14], and 
(ii) if one wishes to monitor the flow of the implementation of an algorithm at any intermediate 
step, then the best option is to measure the full interemediate density matrix [11]. 

For an n-qubit ensemble system, the size of the Hilbert space increases as 2" and the number of 
density matrix elements increases as 2"x2". Of these there are M = (2"~^)(2"~^ -|- 1) independent 
elements, of which (n2"~^) elements are one qubit single quantum observable coherences. To measure 
the remaining elements, a series of one- dimensional experiments with readout pulses to rotate the 
unobservables into observables, have been used [10,11,13,14]. Here we propose a new method for 



quantum state tomography based on the two-dimensional Fourier transform technique, where all 
the off-diagonal elements of a density matrix, both unobservable and observable, are measured in 
a two dimensional experiment. All the diagonal elements are measured in another one-dimensional 
experiment. 

It has been pointed out [25] that the earlier method of tomography [10] involving a large number 
of different measurements; while works well for small spin systems, becomes "prohibitively complex" 
for large spin systems. Such large systems can be easily tomographcd using the proposed method. It 
may be mentioned here, that the proposed method uses a two dimensional experiment and requires 
several ti increments. In principle the size of two dimensional data is independent of the number 
of spins (qubits). However, the number of quantus increase linearly with the number of spins. To 
maintain the same resolution for large number of spins (qubits), the number of ti increments may 
have to be increased appropriately (at best linearly). It may be recalled that the same principle is 
applied to two dimensional NMR of biomolecules, where the size of data is independent of the size 
of biomolecules [27]. The proposed method can also be used for tomography in a wider range of spin 
systems, i.e. quadrupolar or strongly coupled systems, and it requires non-selective r.f. pulses which 
are devoid of errors caused by the selective pulses used by the earher method [10]. The method is 
explained using a 2-qubit system and demonstrated on 2 and 4-qubit systems using simulations. 

II. THE METHOD 

The method is based on the technique of indirect detection of multiple quantum coherences in 
NMR spectroscopy by two dimensional Fourier transform technique [26] , wherein all the off-diagonal 
elements are measured in a two-dimensional experiment (pulse sequence 1(^)). The diagonal ele- 
ments of the density matrix are measured in another one-dimensional experiment (pulse sequence 
1(B)). The NMR pulse sequences for the two experiments are, 

(^4) ti — {TT/'2)y — Gz — Oi-y — detect{t2)i 

(B) G,- Py-detect{t2), (1) 

where ti and t2 are the variable time periods of system evolution, {7:/2)y, a-y and Py are the rf 
pulses, and G^ is the field-gradient pulse. 

In experiment 1(A), a given density matrix (7(0) is allowed to evolve for a time ii, at the end 
of which a (7r/2)j/ pulse transform every element into all other elements of the density matrix, 
including diagonal elements. The G^ pulse dephases the off-diagonal elements averaging them to 



zero, and retains only the diagonal elements. A a^y pulse transforms the diagonal elements into 
all elements of the density matrix including single qubit single quantum coherences. These single 
quantum coherences are then detected as a function of time variable ^2- A series of experiments 
are performed by systematic increment of the ti period and the collected two-dimensional time 
domain data set 5(^1,^2) is double Fourier transformed yielding a two-dimensional frequency domain 
spectrum S'(r2i,r22)- 5'(r2i,f22) contains along 1^2 ah single qubit single quantum coherences and 
along contribution of every off-diagonal elements of the density matrix to these transitions, 
dispersed and displayed by their specific frequency of evolution in the time-domain ti. Cross-sections 
parallel to Qi at one single quantum resonance frequency can be fitted to (7(0), yielding all the off- 
diagonal elements in single two-dimensional experiment. The diagonal elements of cr(0) do not 
contribute to the spectrum. 

To obtain the diagonal elements; experiment 1(B) begins by destroying all off-diagonal elements 
of ct(0) by a gradient pulse, then using a small angle detection pulse to convert difference in diagonal 
elements into observable one qubit single quantum coherences by linear response. The amplitudes 
of the coherences allow calculation of all the off-diagonal elements. The above protocol is explained 
in the following by explicit calculations on a two qubit system. 

III. TWO QUBIT SYSTEM 

Consider a two qubit system consisting of two spin 1/2 nuclei of Larmor frequencies uji and a;2, 
coupled by a weak indirect coupling J. The Hamiltonian for the system is, Ti = uJihz+^^2hz + Jhzhz, 
where Ijz (j = 1, 2) are the spin operators. A selective rf pulse of angle 6 and phase on-resonance 
on spin j corresponds to a unitary transform, f/g^ = exp{—i0(yljxsin(j) -\- IjyCoscj))}. All quantum 
algorithms are implemented in NMR by a specific pulse sequence involving the evolutions under 
system Hamiltonian and rf pulses. A general Hermitian complex trace-less deviation density matrix 
for 2 qubits, has 15 independent elements, is spanned by 15 product operators [27], and can be 
expressed as; 

^Y.^ki hkhu (2) 

k,l 

where k and I can take values 0,1,2,3 corresponding to o, x, y, z respectively, but not simultaneously 
0, and Iio — I20 — I is unit matrix, and q^i are real coefficients. All the elements of the density matrix 
cr(0) can be classified into two groups: (i) diagonal elements involving the deviation populations. 
APfe (deviations from an average population) of various eigenstates and (ii) off-diagonal elements 



involving one-qubit coherences (IQ elements) and multi-qubit coherences (zero and double quantum 
elements for a 2 qubit system). 
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. The trace condition J2k APfc=0 yields 15 independent (12 off-diagonal and 3 diagonal) elements in 
this case. Each IQ element, known as a single quantum element, corresponds to product operators 
of the type Iix, hy, hxhz, hyhz, and I2X, hy, hxhz, and hyhz- Each 2Q element, known as zero 
or double quantum element depending on the frequency of evolution cui —uji or CU1+UJ2, corresponds 
to an expansion in terms of hxhx, hxhy, hyhx and hyhy 



A. Measurement of ofF-diagonal elements 

Effect of the pulse sequence 1(^4) on any initial density matrix a{0) can be described as [27], 
(7(0) ai{h) = e-'^*V(0)e'^*i 

Hmeasure) ^^^^^^ ^ e'^^'^ a,{t,)e^^'\ (4) 

Here the operator V projects and retains only the diagonal part of (T^iti). The complex time domain 
signal s(ii, t-^) obtained on measurement as a function of ^2, is s(ii, — Trace[{Ii ^it)- (75(^1, ^2)], 
which after double Fourier transform gives two-dimensional (2D) spectrum S{yLi,VL2) which is a 
function of the two frequency variables Vti and VL2. The ^2 axis of this spectrum has only the 
single quantum (IQ) elements (four transitions in the 2-qubit system). Along the VLi axis of the 2D 
spectrum all the off-diagonal elements of (7(0) yield peaks corresponding to their amplitudes (g^;) 
and individual evolution frequencies during the time period ti. 

Cross-sections of the signal after Fourier Transform with respect to t2, but before ti, [S{ti,Q,i)\ 
taken parallel to ti at ^2 = oj^ and ^2 = <^2^ respectively correspond to, 

ai{u^) = e'^''/'^'^sina[{qiQ qi^)cos{utti) {q^ - qiz)cos{u^ti) 



-(^20 + q23)sin{ujfti) - (§20 - fejs^l'^r^i) 
+^cosa{{qu - q22)cos{u^ti) + {qu + q22)cos{uj^ti) 

(gi2 + q2i)sin{uj^ti) - {qu - q2i)sin{uj^ti)}]/4: 



(74(0;^) — e 



■tl/T2 



sina[{qQi + q2.i)cos{ujtti) + (goi - q3i)cos{uj^ti) 
-{qo2 + q32)sin{uj^ti) - (go2 - q32)sin{u2ti) 

+^cosa{{qu - q22)cos{uj^ti) + {qu + ^22 ) cos (tu^ti) 
{qi2 + q2i)sin{uj^ti) - {q^ - q2i)sin{uj^ti)}]/4:, 



(5) 



where uj 



Uj ± J/2, uj- 



uii + oj2i ui^ = coi — UJ2, qki are the coefficients of expansion of a{0) as 



in Eq. (2), and T2 is the transverse relaxation time of various coherences. T2 can be different for 
each coherence, but taken identical here for simplicity. On ti Fourier transformation of expression 5, 
one obtains the two-dimensional (2D) spectrum S{fli, O2) in which the cosine terms give absorptive 
and the sine terms dispersive Lorenztian lines. All the coefficients of off-diagonal elements of a{0) 
can be obtained by fitting the cross-sections from the two-dimensional spectrum S{D,i,D,2) (taken 
parallel to Qi) to the absorptive/ dispersive Lorentzians obtained from expressions similar to Eq. (5). 
Two cross-sections, one at each qubit, are sufficient to calculate all the off-diagonal elements of 
the 2-qubit density matrix cr(0) of Eq. (2). For an n-qubit system, there are 2"~^ cross-sections 
per qubit. While only one cross-section per qubit is required to map all the off-diagonal elements, 
the remaining cross-sections can be used to minimise the errors. It may also be noted that the 
diagonal elements do not interfere with the two-dimensional spectrum obtained by pulse sequence 
1(A). Similarly the undesired 'axial peaks' (having zero frequency during ti period) arising due to 
the longitudinal relaxations during ti period [27] are also suppressed by the present scheme. 

In the two-qubit case, the one-qubit and two-qubit coherences have different conversion ratios, 
namely, sina and \sin2a. As a compromise a — tt/A has been used here. For qubit systems having 
higher-qubit coherences, appropriate values of a should be used which optimizes the intensities of 
various orders [28]. 



B. Measurement of diagonal elements 



The transformations of the density matrix a{0) by the pulse sequence 1{B) are as follows. 



a{0) 



q3ohz + qozhz + qmhzhz 




qm{hzCos(3 IixSinP) + qQz{l2zCos(5 + hxSinP) 



+q33{hzCosP + hxSin(3){l2zCos(3 + hxSinP) 

Measure. (6) 

To measure the coefficients under linear response, /? should be small [27]. The measured one-qubit 
coherences then are 

PiQsohx + Qoshx + ?33 {hxhz + hzhx)] , (7) 

which can be rearranged as, 

/3[(g30 + q33/2)ihx + 2hj2z) + (gao - q^^/'^Wix - ^hxhz) 

+ {Q03 + q33/2){l2x + 2hzhx) + (g03 " q33/2){l2x " 2hj2x)]/'2. (8) 

The coefficients of the four terms in the above expression are proportional to the intensities of the 
corresponding four transitions of a two qubit system. After calculating g'30, qos, and g'33, all the 
diagonal elements of the density matrix a{0) can be calculated, since the diagonal part is equal to 

qsohz + qoshz + qsahzhz- 

It should be noted that the gradient pulse Gz used in pulse sequences 1(A) to destroy off-diagonal 
elements, does not destroy homonuclear zero-quantum coherences. In such cases, an extra small 
delay along with G^, randomly varied between each ti experiment can suppress the homonuclear 
zero quantum coherence [27]. In experiment 1(B), signal averaging using a few randomly varied 
along with G^ would suppress the homonuclear zero quantum coherences. 

The above schemes 1(A) and (B) assume ideal r.f. pulses. To correct for errors due to imper- 
fection of the r.f. pulses a third one dimensional experiment can be performed to measure the one 
qubit coherences directly without application of any pulses. These coherences can then be used to 
normalize all other elements of the density matrix measured by experiments 1(A) and (B). 

IV. SIMULATION 

To demonstrate the protocol we tomograph an arbitrary complex density matrix with simulations. 
In a 2-qubit system, such a density matrix is of the form; 

(t{0) = hz + 2.3/2. + 6.77^/2. + hx + lOhxhz + 5hy + S.Bhyhz + 2.5Ilyl2y 

+ 7.2hyhx + ^^hxhx + 1.45/1^/22, + 2/20. + ^Abhzhx + 6.9/23, + Q-75^hzhy 



/ 3.325 1.8625 - 5. 1383i 3 - 3.75i 2.65 - 2.1625?: \ 

1.8625 + 5.1383i -2.325 3.875 - 1.4375i -2 - 1.625i 

(9) 

3 + 3.375i 3.875 + 1.4375i -1.025 0.1375 - 1.7618i 

V 2.65 + 2.1625i -2 + 1.625i 0.1375 + 1.7618i 0.025 / 

We assume the Larmor frequencies of the two qubits (spins) as a;i=1200 Hz, a;2=1800 Hz and the 
indirect coupUng constant between the qubits as J=200 Hz. Experiment 1(A) is performed to obtain 
all the off-diagonal elements, a was chosen as 45". 512 ti increments were performed yielding the 
2D spectrum shown in Fig. 1(a). Cross sections parallel to taken at one of the transitions of 
each qubit are shown in Fig 1(b) and (c). These cross sections were fitted to get all the complex 
off-diagonal elements of the density matrix. T2 was taken as 10ms for all coherences. The diagonal 
elements were mapped using experiment 1(B) with /3 — 10°, Fig 1(d). The real and imaginary parts 
of the tomographed density matrix arc shown respectively in Figs. 1(c) and 1(f). The calculated 
density matrix matches the input density matrix better than 0.01% for all complex elements. We 
have also carried out the simulations on a 4-qubit system (Fig. 2) and tomographed the density 
matrix with 99.7% fidelity. 



V. CONCLUSION 



Two-dimensional nuclear magnetic resonance spectroscopy provides an efficient method for the 
quantum state tomography. Only an one-dimensional and a two-dimensional experiment are required 
for measuring all the elements of the density matrix. Since the earlier method requires a series of one 
dimensional experiments with different readout pulses, for large spin systems the approach becomes 
enormously complex [16]. However such systems can easily be tomographed using the proposed 
method by aptly increasing the ti increments. The earlier method uses spin-selective r.f. pulses, 
which requires long-duration irradiation of a paricular spin. During such a pulse, the unperturbed 
spins evolve under the Zeeman and coupling interactions, introducing errors due to measurement 
[15,16]. The method described here requires non-selective short-duration r.f. pulses which do not 
introduce such errors. Search of more qubits has led researchers to use strongly coupled spin-1/2 
nuclei and quadrupolar nuclei (spin>l/2) oriented in liquid crystalline matrices [19-23]. For such 
systems the notion of spin-selectivity does not apply, but the proposed method based on non-selective 
pulses can be used for tomography. Recently, the method was used to tomograph the states while 
quantum informtion processing in weakly and strongly coupled spin systems [23,24]. This method 
can also be extended to a 3-dimensional experiment in which quantas of various orders are displaced 



in different planes of the 3-D experiment, increasing the detectibihty and the resolution of the 
spectrum [29]. 
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FIG. 1. Tomography of a complex density matrix in a 2 qubit system, (a) is the 2D spectrum generated by experiment 
1(A) for tomography of all off diagonal elements, (b) is the cross section parallel to Sli taken at Q2= 1300 Hz {uji + J/2) 
(a transition frequency of first qubit), (c) is the cross section parallel to Sli taken at i^2— 1900 Hz { UJ2 + J/'2) (a transition 
frequency of the second qubit) and (d) is the ID spectrum obtained by experiment of 1(B) for mapping all diagonal elements of 
(t(0). (e) and (f) are the real and imaginary part of the tomographed density matrix. The calculated values of the tomographed 
matrix are an = 3.2501, aw = 1.8625 - 5.1385i, ars = 3.0001 - 3.375i, ai4 = 2.6501 - 2.1624i, a2i = 1.8625 + 5.1385i, 
a22 = -2.3251, (723 = 3.8750 + 1.4374i, a24 = -2.0001 - 1.625i, (J31 = 3.0001 + 3.375j, (732 = 3.8750 - 1.4374i, (J33 = -1.025, 
C734 = 0.1376 - 1.76 18i, a^i = 2.6501 + 2.1624i, a42 = -2.0001 + 1.625i, a43 = 0.1376 + 1.7618i, and (744 = 0.025. 



FIG 2. Tomography of a complex density matrix in 4-qubit system. An arbitrary complex density matrix cr(0) — 0.8/ia; + Iiy + 

0. 512X + l2y + 0.9/32; + 1.1/31; + ^4a; + 1.2/4y + 6.3Ilx l2x hxlix + S.9Ilx l2y Isy l4y + Ilxl2xl3zl4z + ^■H2xl'ixlix + ^■Qllxl2xlzy + 

1. ^I\xl2zHyIix + 0.6/42 + l2zliz + 1.3/2z/3z/4z + '2Iizl2zI-izI'iz , is reconstructed with the 2D Fourier Transform technique. The 
frequencies and couplings used in the simulation are uji — 600Hz, UJ2 = 750Hz,ujs = 1000Hz,uj4 — UOOHz, J12 — 20// 2:, J13 = 
10Hz,Ji4 = 7O//2, J23 = 35//z, J24 = 2AHz and J34 = 16Hz. (a) Shows the real part of (7(0), (b) is the real part of 
reconstructed density matrix and (c) the differences ( magnified by 10"^) between the elements of (a) and (b). The 256 complex 
elements of (t{0) (imaginary part not shown here) were tomographed with more than 99.7% accuracy. 



